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Abstract 



We consider a self-gravitating ideal fermion gas at nonzero tem- 
perature as a model for the Galactic halo. The Galactic halo of mass 
O I ~ 2 X IO^^Mq enclosed within a radius of ~ 200 kpc is consistent with 

the existence of a supermassive compact dark object at the Galactic 
■ center that is in hydrostatic and quasi-thermal equilibrium with the 

halo. The central object has a maximal mass of ~ 2.3 x 10^ within 
^ I a minimal radius of ~ 18 mpc for fermion masses ~ 15 keV. 

H , 

. >r , In the past, self-gravitating neutrino matter was suggested as a model for 

quasars, with neutrino masses in the range 0.2keV < m < O.SMeV More 
recently, supermassive compact objects consisting of nearly non-interacting 
degenerate fermionic matter, with fermion masses in the range 10 < m/keV 
^ 20, have been proposed p|, |^, ||, ^, |^ as an alternative to the supermassive 
black holes that are believed to reside at the centers of many galaxies. 

So far the masses of ~ 20 supermassive compact dark objects at the 
galactic centers have been measured 0. The most massive compact dark 
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object ever observed is located at the center of M87 in the Virgo cluster, and 
it has a mass of ~ 3 x lO^M© If we identify this object of maximal mass 
with a degenerate fermion star at the Oppenheimer-Volkoff (OV) limit [0], 



Mov = 0.54M|im-2^-i/2 ^ 3 ^ lO^M© g, where Mpi = ^hc/G, this 
allows us to fix the fermion mass to m ~ 15 keV for a spin and particle- 
antiparticle degeneracy factor g = 2. Such a relativistic object would have a 
radius i?ov = 4.45i?s — 1-5 light-days, where -Rs is the Schwarzschild radius 
of the mass Mqv- It would thus be virtually indistinguishable from a black 
hole of the same mass, as the closest stable orbit around a black hole has a 
radius of 3 -Rs anyway. 

Near the lower end of the observed mass range is the compact dark object 
located at the Galactic center with a mass — 2.6 x IO^Mq. Interpret- 
ing this object as a degenerate fermion star consisting of m ~ 15 keV and 
g = 2 fermions, the radius is Rc — 21 light-days ~ 7 x 10^i?s 0, -Rs being 
the Schwarzschild radius of the mass Mc- Such a nonrelativistic object is far 
from being a black hole. The observed motion of stars within a projected 
distance of ~ 6 to ~ 50 light-days from Sgr A* [jlOl yields, apart from the 
mass, an upper limit for the radius of the fermion star -Rc ^ 22 light-days. 

The required nearly non-interacting fermion of ~ 15 keV mass cannot be 
an active neutrino, as it would overdose the Universe by orders of magnitude 
[ p!T[] . However, the ~ 15 keV fermion could very well be a sterile neutrino 
mixed to active neutrinos with a mixing angle sin^ 2^ ~ 10~^^. Indeed, as 
has been shown for an initial lepton asymmetry of ~ 10~^, a sterile neutrino 
of mass ~ 10 keV may be resonantly or non-resonantly produced in the early 
Universe with near closure density, i.e., Qd ~ 0.3 [Q. As an alternative 
possibility, the required ~ 15 keV fermion could be the axino [|l^ or the 
gravitino in soft supersymmetry breaking scenarios. 

In the recent past, galactic halos were successfully modeled as a self- 
gravitating isothermal gas of particles of arbitrary mass, the density of which 
scales asymptotically as r~^, yielding flat rotation curves ||15|. As the super- 
massive compact dark objects at the galactic centers are well described by 
a gas of fermions of mass m ~ 15 keV at T = 0, it is tempting to explore 
the possibility that one could describe both the supermassive compact dark 
objects and their galactic halos in a unified way in terms of a fermion gas at 
finite temperature. We show in this letter that this is indeed the case, and 
that the observed dark matter distribution in the Galactic halo is consistent 
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with the existence of a supermassive compact dark object at the center of 
the Galaxy which has about the right mass and size. 

Degenerate fermion stars are well understood in terms of the Thomas- 
Fermi theory applied to self-gravitating fermionic matter at T = p|. Ex- 



tending this theory to nonzero temperature [0, |T^, , it has been shown 
that at some critical temperature T = T^, a. self-gravitating ideal fermion 
gas, having a mass below the OV limit enclosed in a sphere of radius R, may 
undergo a first-order gravitational phase transition from a diffuse state to 
a condensed state. However, this first-order phase transition can take place 
only if the Fermi gas is able to get rid of the large latent heat. As short-range 
interactions of the fermions are negligible, the gas cannot release its latent 
heat; it will thus be trapped for temperatures T < Tc in a thermodynamic 
quasistable supercooled state close to the point of gravothermal collapse. 
The formation of a supercooled state close to the point of gravothermal 



collapse may be due to violent relaxation |]20|, Through the gravi- 

tational collapse of an overdense fluctuation, ~ 1 Gyr after the Big Bang, 
part of gravitational energy transforms into the kinetic energy of random 
motion of small-scale density fluctuations. The resulting virialized cloud 
will thus be well approximated by a gravitationally stable quasi-thermalized 
halo. In order to estimate the particle mass-temperature ratio, we assume 
that a cold overdense cloud of the mass of the Galaxy M, stops expanding 
at the time t^, reaching its maximal radius and minimal average density 
Pm = 3M/(47ri?^). The total energy per particle is just the gravitational 
energy 

E=-'-^. (1) 

5 R^a 



Assuming spherical collapse p3|, one arrives at 



9n^ _ 9it'^ 

Pm = —P{tm) = — ^^dPo(l + ZraY, (2) 

io ib 

where p(tm) is the background density at the time or the cosmological 
redshift z^, and po = 3Hq/{StiG) is the present critical density. We now 
approximate the virialized cloud by a singular isothermal sphere of mass 
M and radius i?, characterized by a constant circular velocity G = (2T/m)^/^ 
and the density profile p(r) = 6^/47rGr^. Its total energy per particle is the 
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sum of gravitational and thermal energies, i.e., 

IGM 1 2 



Combining Eqs. (P, (|[), and (|]), we find 

e' = ^G{6QdPoM'f'{l + z^). (4) 

Taking = 0.3, M = 2 x lO^^^vf^^ ^ 4^ ^nd Hq = 65kms-^Mpc-\ 
we find G ~ 220kms~^, which corresponds to the mass-temperature ratio 
m/T ~ 4 X 10^. 



We now briefiy discuss the Thomas-Fermi theory [T^, |1^ for a self- 
gravitating gas of N fermions with mass m at the temperature T enclosed 
in a sphere of radius R. We restrict our attention to the Newtonian theory 
since the general relativistic effects are not relevant to the physics we discuss 
in this paper. The general relativistic treatment will be reported elsewhere. 
For large A^, we can assume that fermions move in a spherically symmetric 
mean-field potential ip{r) which satisfies Poisson's equation 

^ = G^; !^ = 4.rW. (5) 

dr dr 

Ai being the enclosed mass. The number density of fermions (including 
antifermions) n can be expressed in terms of the Fermi-Dirac distribution (in 
units h = c = k = 1) 



p f d^q ( q 



2 



m 



-1 



Here g denotes the combined spin-degeneracy factor of the neutral fermions 
and antifermions, i.e., is 2 or 4 for Majorana or Dirac fermions, respectively. 
For each solution ip{r) of (||), the chemical potential p is adjusted so that the 
constraint 

R 

drA'nr'^n{r) = iV, (7) 



is satisfied. Equations with (H) should be integrated using the boundary 
conditions at the origin, i.e., 

<^(0) = <^o; A^(0) = 0. (8) 
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It is useful to introduce the degeneracy parameter 



V = j;-Y^- (9) 

As is monotonously increasing with increasing r, the strongest degeneracy 
is obtained at the center with r/o = {fJ'~^Vo)/T. The parameter rjo, uniquely 
related to the central density, will eventually be fixed by the constraint (|^) 
or equivalently by the condition A4{R) = mN at the outer boundary. In 
this way, the explicit dependence on the chemical potential fi is absorbed in 
the degeneracy parameter r]o. For r > R, the function (f yields the usual 
empty-space Newtonian potential 

^{r) = -—-. (10) 

The set of self-consistency equations with the boundary conditions 

(H), defines the gravitational Thomas- Fermi equation. 

The numerical procedure is now straightforward. For a fixed, arbitrarily 
chosen ratio m/T, we first integrate Eqs. (§) numerically on the interval 
[0, i?] to find the solutions for various central values ?7o. This yields M.{K) 
as a function of ?7o. We then select the value of ?7o for which }A[K) = mN. 

The quantities A^, T,and R are free parameters in our model and their 
values are dictated by physics. In the following, N is required to be of 
the order 2 x lO^^M©/??!, so that for any m, the total mass is close to the 



estimated mass of the halo As we have demonstrated, the expected 

temperature of the halo is given by m/T = 4 x 10^. Our choice R = 200 kpc 
is based on the estimated size of the Galactic halo. The only remaining free 
parameter is the fermion mass, which we fix at m = 15 keV, and justify its 
choice a posteriori. 

For fixed A^, there is a range of T where the Thomas-Fermi equation has 
multiple solutions. For example, for N = 2 x 10^^ and m/T = 4 x 10^, 
we find six solutions, which we denote by (a), (b), (c), (c'), (b'), and (a') 
corresponding to the values rjo = -30.53, -25.35, -22.39, 29.28, 33.38, and 
40.48, respectively. In Fig. |l| we plot the mass density profiles of the halo. 
For the negative central value rjo, for which the degeneracy parameter is 
negative everywhere, the system behaves basically as a Maxwell-Boltzmann 
isothermal sphere. Positive values of the central degeneracy parameter r/o 
are characterized by a pronounced central core of mass of about 2.5 x lO^M© 
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within a radius of about 20 mpc. The presence of the core is obviously due 
to the degeneracy pressure. The core represents material which, having been 
cooled by expansion, experiences little entropy increase during the ensuing 
collapse. Thus the dynamics of its formation should be well approximated 
by a dynamical Thomas-Fermi theory based on the equation of state of a 
degenerate Fermi gas [^. Conversely, the halo is formed from phase- mixed 



matter and estimates similar to those leading to (^) give an average entropy 
per particle increasing from few x 10" to few x 10^. 

A similar structure was obtained in coUisionless stellar systems modeled 
as a nonrelativistic Fermi gas [B^. Note that while violent relaxation leads 



to a Fermi-Dirac distribution in either case, for stars the onset of degeneracy 
signals the breakdown of the assumption that collisions are unimportant, re- 



sulting in a Maxwell-Boltzmann distribution . No such breakdown occurs 
for elementary fermions [^. 

Fig. |I| shows two important features. First, a galactic halo at a given 
temperature may or may not have a central core depending on whether the 
central degeneracy parameter rjo is positive or negative. Second, the closer 
to zero ?7o is, the smaller the radius at which the asymptotic behavior 
of density begins. The flattening of the Galactic rotation curve begins in 
the range 1 ^ r/kpc ^ 10, hence the solution (c') most likely describes the 
Galactic halo. This may be verified by calculating the rotation curves in 
our model. We know already from our estimate (^) that our model yields 
the correct asymptotic circular velocity of 220 km/s. In order to make a 
more realistic comparison with the observed Galactic rotation curve, we must 
include two additional matter components: the bulge and the disk. The bulge 



is modeled as a spherically symmetric matter distribution of the form [2S 



Q-hs »oo ^-hsu 

^"'^'^ = ^1 ^^ [(^ + i)«-i]V2 ' (11) 

where s = (r/roY^^, tq is the effective radius of the bulge and h is a parame- 
ter. We adopt tq = 2.67 kpc and h yielding a bulge mass Mb = 1.5 x 10^°Mq 
. In Fig. 1^ the mass of the halo and bulge enclosed within a given radius 
is plotted for various rjo. The data points, indicated by squares, are the mass 
Mc = 2.6 X IO^Mq within 18 mpc, estimated from the motion of the stars 
near Sgr A* ]l^, and the mass M50 = 5.4l3g x 10^"^ within 50 kpc, estimated 



from the motion of satellite galaxies and globular clusters p5[. Variation of 
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Figure 1: The mass density profile of tlie lialo for ?7o = (dotted line) and 
for the six r^Q- values discussed in the text. Configurations with negative rjQ 
((a), (b), (c)) are depicted by the dashed and those with positive rjQ ((a'), 
(b'), (c')) by the solid line. 
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Figure 2: Enclosed mass of the halo plus bulge versus radius for rjo = 24 
(dashed), 28 (solid), and 32 (dot-dashed line). 
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the central degeneracy parameter rjQ between 24 and 32 does not change the 
essential halo features. 

In Fig. 1^ we plot the circular velocity components: the halo, the bulge, 



and the disk. The contribution of the disk is modeled as 30 



where we take Tq = 13.5 kpc and Gd = 100 km/s. Here we have assumed for 
simplicity that the disc does not influence the mass distribution of the bulge 
and the halo. Choosing the central degeneracy r/o = 28 for the halo, the data 



by Merrifield and Oiling |31| are reasonably well fitted. 

We now turn to the discussion of our choice of the fermion mass m = 15 
keV for the degeneracy factor g = 2. To that end, we investigate how the 
mass of the central object, i.e., the mass within 18 mpc, depends on m 
in the interval 5 to 25 keV, for various tjq. We find that m ~ 15 keV always 
gives the maximal value of ranging between 1.7 and 2.3 xIO^Mq for f]Q 
between 20 and 28. Hence, with m ~ 15 keV we get the value closest to the 
mass of the central object Mc estimated from the motion of the stars lear 



Sgr A* P 



The radius of our central object of about 18 mpc is much larger than the 
size of the radio source Sgr A* . In fact , very large array interferometric obser- 
vations of Sgr A* at millimeter wavelength show that the radiowave emitting 



region is < 1-3 AU [^. However, it has not yet been shown conclusively 
that Sgr A* is indeed the object that has a mass ~ 3 x 10^ Mq. There are 
arguments, based on the nonmotion of Sgr A* and equipartition of energy 
in the central star cluster, indicating that Sgr A* could have a mass > 10^ 
Mq . This argument is only conclusive if equipartition of energy actually 
takes place in a reasonable time frame. It is, therefore, still possible that the 
compact radiosource Sgr A*, with a radius of a few AU, and the moderately 
compact supermassive dark object that has been detected gravitationally, 
and possibly also in X-rays in the quiescent state, with a radius of ~ 20 mpc 
[0, could be two distinct objects. 

In summary, using the Thomas-Fermi theory, we have shown that a 
weakly interacting self-gravitating fermionic gas at finite temperature yields 
a mass distribution that successfully describes both the center and the halo 
of the Galaxy. For a fermion mass m ^ 15 keV, a reasonable fit to the 
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Figure 3: Fit to the Galactic rotation curve. The data points are by Olhng 



and Merrifield for = 8.5 kpc and Gq = 220 km/s. 
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rotation curve is achieved with the temperature T = 3.75 meV and the 
degeneracy parameter at the center r]o = 28. With the same parameters, 
the masses enclosed within 50 and 200 kpc are M50 = 5.04 x IO^^Mq and 
^200 = 2.04 X IO^^Mq, respectively. These values agree quite well with the 
mass estimates based on the motion of satellite galaxies and globular clus- 
ters [^. Moreover, the mass Mc ~ 2.27 X IO^Mq, enclosed within 18 mpc. 



agrees reasonably well with the observations of the compact dark object at 
the center of the Galaxy. We thus conclude that both the Galactic halo and 
the center could be made of the same fermions. 
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